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An algorithm is presented which determines the equivalence of two settings of a 
(3 + d) -dimensional superspace group (d = 1, 2, 3). The algorithm has been 
implemented as a web tool finds sg on SSG(3 + d)D, providing the 
transformation of any user-given superspace group to the standard setting of 
this superspace group in SSG(3 + d)D. It is shown how the standard setting of a 
superspace group can be directly obtained by an appropriate transformation of 
the external-space lattice vectors (the basic structure unit cell) and a 
transformation of the internal-space lattice vectors (new modulation wavevec- 
tors are linear combinations of old modulation wavevectors plus a three- 
dimensional reciprocal-lattice vector). The need for non-standard settings in 
some cases and the desirability of employing standard settings of superspace 
groups in other cases are illustrated by an analysis of the symmetries of a series 
of compounds, comparing published and standard settings and the transforma- 
tions between them. A compilation is provided of standard settings of 
compounds with two- and three-dimensional modulations. The problem of 
settings of superspace groups is discussed for incommensurate composite 
crystals and for chiral superspace groups. 



1. Introduction 

Symmetry is one of the most important concepts in the solid- 
state sciences. Knowledge of the symmetry of a crystalline 
compound allows the understanding of many aspects of its 
physical behavior, including degeneracies, the possibility of 
possessing non-linear properties and the anisotropy of the 
response to external fields. A change in symmetry at different 
temperatures, pressures or compositions is used as the key 
parameter for characterizing phase transitions of a compound. 
Symmetry is used for the description of phonon and electron 
bands and thus allows the interpretation of spectroscopic 
measurements on materials. Not least, symmetry restrictions 
on structural parameters are essential for successful refine- 
ments of crystal structures. 

Theoretically, the classification of symmetry is solved. The 
230 space groups give the 230 possibilities for the symmetry of 
a periodic structure (Hahn, 2002). Aperiodic crystals lack 
three-dimensional (3D) translational symmetry (Janssen et al, 
2006, 2007; van Smaalen, 2007). The structures of incom- 
mensurately modulated crystals are characterized by a three- 
dimensional lattice for the average structure together with d 
modulation waves (d = 1,2, . . .) describing deviations from 
the lattice-periodic structure. Their symmetries are given by 
(3 + d)-dimensional [(3 + d)Y)] superspace groups (de Wolff et 
al, 1981). The latter are space groups of (3 + d)D space, which 
have to obey particular conditions in order to qualify as 



symmetry groups for the symmetries of aperiodic crystals. 
Incommensurate composite crystals are described by the same 
superspace groups as modulated crystals (van Smaalen, 2007), 
while quasicrystals require a slightly modified treatment 
(Janssen et al, 2007; Steurer & Deloudi, 2009). 

We have recently generated a complete list of superspace 
groups and their Bravais classes of dimensions (3 + d) for d = 
1, 2 and 3 (Stokes et al, 2011a). The list agrees with previous 
information on (3 + 1)D superspace groups (Janssen et al. , 
2006), but it contains numerous corrections for superspace 
groups of dimensions d = 2 and 3 (Yamamoto, 2005) and even 
some corrections to the Bravais classes of dimensions 2 and 3 
(Janssen et al, 2006). The results of Stokes et al (2011a) are 
compiled in the form of the web-based data repository 
SSG(3 + d)D (Stokes et al, 20116). SSG(3 + d)D provides 
several types of information for each superspace group, 
including the Bravais class, the list of symmetry operators and 
reflection conditions in both standard and supercentered 
settings. 

It is noticed that Stokes et al (2011a) have defined the 
standard settings and their symbols by a set of judiciously 
chosen rules, which, however, include subjective choices. The 
standard setting thus is defined as the setting included in the 
list of superspace groups on SSG(3 + d)D. 

The use of alternate settings of space groups is a well 
known feature for three-dimensional space groups. Volume 
A of the International Tables for Crystallography (Hahn, 2002) 
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provides several settings for monoclinic space groups, thus 
showing, for example, that C2/c, A2/n and 12/ a denote 
different settings of space group No. 15. Alternate settings of 
three-dimensional space groups arise owing to different 
choices of the unit cell. Trivial transformations include a 
simple relabeling of the axes a, b and c. For monoclinic space 
groups this implies the freedom to select the unique axis as 
a = a l5 b = a 2 or c = a 3 . For the example of space group No. 
15 the transformation 

a l = a l + a 3> 3 2 = a 2> a 3 = — a l (1) 

takes the setting C2/c into A2/n, while both settings refer to a 
unique axis a 2 . Applying the transformation of equation (1) to 
the setting A2/n results in the third setting 12/ a. Other 
notorious pairs of equivalent settings include the / and F 
settings for space groups based on the /-centered tetragonal 
lattice, primitive and centered hexagonal settings as well as 
obverse versus reverse settings for rhombohedral space 
groups, and the //-centered setting as an alternative to the 
primitive setting for trigonal space groups (Hahn, 2002). 

Different diffraction experiments independently lead to any 
of the possible settings of the space group. It is then of high 
practical importance to find the transformation between these 
settings or to establish that two different space groups have 
indeed been obtained. The latter situation implies different 
compounds or different phases of one compound, while 
different settings of one space group imply that the same 
compound has been studied. In other experiments it is 
important that a previously defined orientation of a crystalline 
material is re-established, thus requiring the relation to be 
found between the newly found setting and a standard setting. 

As a start, superspace groups for incommensurately mod- 
ulated compounds and incommensurate composite crystals 
exhibit the same variation of settings as three-dimensional 
space groups do, since (3 + d)D superspace groups are based 
on a basic structure lattice and space group in three- 
dimensional space. In addition, superspace groups may appear 
in many more different settings, owing to the ambiguity in the 
choice of the modulation wavevectors characterizing the 
structure and the diffraction pattern. The equivalence of 
different settings of a superspace group is not always obvious. 
In some cases, establishing an actual equivalence can be a 
computationally prohibitive task unless appropriate algo- 
rithms are used. Here we present such an algorithm, which was 
used but not described in detail in our previous publication 
Stokes et al. (2011a). It is available within SSG(3 + d)D as a tool 
with which to determine the transformation between a user- 
provided setting of a (3 + d)D superspace group and the 
standard setting defined by SSG(3 + d)D (Stokes et al, 2011b). 
It thus can be used to establish or disprove the equivalence of 
settings. 

Coordinate transformations between different settings of 
three-dimensional space groups are discussed in Volume A of 
International Tables for Crystallography (Hahn, 2002). For 

1 The mathematically correct designation is that of equivalent space groups 
providing different settings of space-group type No. 15. 



(3 + 1)D superspace groups, typical transformations are 
presented in van Smaalen (2007). The possibility to combine 
two modulation wavevectors into an equivalent but different 
set of two wavevectors leads to new types of transformations 
for d = 2 and 3. 

One goal of this paper is to present an overview of typical 
coordinate transformations that may occur between settings of 
superspace groups. Particular attention is given to the relation 
between the formal description in (3 + d)D space, as given on 
SSG(3 + d)D (Stokes et al, 2011a), and an experimentally 
related description in terms of a rotation and an origin shift in 
three dimensions (van Smaalen, 2007). Where available, we 
use substances published in the literature to illustrate impor- 
tant transformation types. 



2. Equivalence of superspace groups 

2.1 . Definitions 

The following definitions are used by Stokes et al. (2011a), 
van Smaalen (2007) and Janssen et al. (1995). A d-dimen- 
sionally modulated structure is characterized by d rationally 
independent modulation wavevectors q 7 with components 
compiled in a d x 3 matrix a according to 

q , ' = o) 1 a 1 +ff ;2 a2+ff ;3 a;. (2) 

For an aperiodic structure at least one component in each 
row of a is an irrational number. The reciprocal vectors 
{af, a|, a 3 , q 1 , . . . , q d ] in physical space correspond to the 
reciprocal basis vectors {a* lt . . . , a* 3+d } in superspace. The 
basis vectors of the direct lattice in superspace are 
{a sl , . . . , a s3+rf } and the coordinates of a point in superspace 
are (x sl , ...,x s3+d ). 

Note that the SSG(3 + d)D data repository, and also the web- 
based finds sg and transforms sg tools described herein, 
presently use the {x,y,z,t\ notation to indicate superspace 
coordinates, though the same notation is also commonly used 
to indicate physical- and internal-space coordinates. 

An operator g of a (3 + d)D superspace group G consists of 
a rotation R s (g) and a translation v s (g) given in matrix form as 

«•>=(; ;). =(::). « 

where R is a 3 x 3 integer matrix and v 3 is a three-dimensional 
column vector, together defining the operator {Z?|v 3 } in 
physical space, e is a d x d integer matrix, and det(Z?) = 
det(s) = ±1. The d x 3 integer matrix M is defined as 
[equation (2)] 

M = oR i - so. (4) 

M has nonzero components only in the case that at least one of 
the modulation wavevectors incorporates nonzero rational 
components. Following Stokes et al. (2011a), each operator g 
can be written as an augmented (4 + d) x (4 + d) matrix 
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A(g) = 




(5) 



that simultaneously treats the point and translational parts of 
the operation. The action of operator g on a point in super- 
space then is given by the matrix product 



V 



- s3+d 
1 



A(g) 



( X sl \ 



\ 1 / 



(6) 



In superspace a coordinate transformation can be accom- 
plished by the augmented affine transformation matrix 



S s S v 
0 1 




\d 



(7) 



The components of S R , S M and S e are required to be integers. 
Also, det(S s ) = 1 and det(5' E ) = ±1. The transformation S can 
be interpreted as a rotation S s in superspace followed by a 
change of origin S v . The effect of this transformation in 
physical space can be described in terms of a rotation of the 
reciprocal basis and the choice of an alternate set of modu- 
lation wavevectors, according to 



a i* — E ^Rik a k 

k=l 

a 1 = {S M + S e o)S- R \ 



(8) 



Two (3 + d)D superspace groups, G 1 and G 2 , are equivalent if 
a single transformation S can be found, such that for every 
g 2 G G 2 , 



Atf) = SA(g 2 )S- 



(9) 



for some g 1 e G 1 . Note that these definitions imply that a 
primitive setting is used for the superspace group, where all 
lattice translations are represented by integers, even those 
which are centering translations in a conventional setting 
(Stokes et ai, 2011a). 

It is sufficient that the relation of equivalence [equation (9)] 
is tested for corresponding pairs of non-translational genera- 
tors from the two superspace groups; the generators of the 
translation subgroup need not be considered. Furthermore, a 
transformation of the type of equation (9) can only be found 
between operators g 11 and g 21 if [equation (3)] 



det(/?") = det(R 2 '), 
det(e") = det(s 2 '), 



trace(i? l; ) = trace(i? 2/ ), 
trace(e 1 ') = trace(e 2 '), 



(10) 



where g is the /th generator of G 1 and g 21 is the corresponding 
generator of G 2 . The appropriate pairs of generators are 
obtained by consideration of the basic structure space group 
implied by the superspace group. Other, trivial properties that 
need to be fulfilled for equivalence and that are easily tested 
include the number of operators in the point group of the 
superspace group, which must be equal for G 1 and G 2 . 



2.2. The algorithm determining equivalence 

The goal of testing for equivalence of two superspace 
groups is to find the augmented matrix S with which the 
operators of G 2 are transformed into corresponding operators 
of G 1 , or to establish that a matrix S that solves equation (9) 
simultaneously for all pairs of generators does not exist. 
Equation (9) is quadratic in S but can be recast in linear form 
as SA(g 2 ) = A(g 1 )S. Given that a pairing has been established 
for the « gen generators of G 1 and G 2 , this results in (3 + dfn, 
equations for (3 + df variables S sik , 

3+d 

E s S ijRsjk ~ S s jk = o, 



gen 



(11) 



and (3 + d)n gen equations for additional (3 + d) variables S vi , 

3+d 3+d 

S vi - E Rsl}S vi = y s f ~ E S sii y s f (mod 1). (12) 

;=i ;=i 

The translational parts of the operators g u and g 21 are only 
known up to a lattice translation, which is taken into account 
by the mod 1 in equation (12). 

Employing the special structure of S [equation (7)], the 
variables S sik can be ordered in a column vector as 



Sr 



(13) 



where, for example, S M is obtained by juxtaposition of the 
columns of S M into a single column matrix. This procedure 
eliminates the 3d variables that are zero according to 
equation (7) and results in [(3 + df — 3d]n gen equations in 
[(3 + df — 3d] variables S sj , 



(3+d) 2 -3d 

E 

7=1 



B tj S sj = 0. 



(14) 



The [(3 + df - 3d]n gen x [(3 + df - 3d] matrix B tj is 
obtained by rearrangement of equation (11), followed by 
linear row operations that bring it into row echelon form. In 
this form, the first nonzero element in each row occurs in a 
column where it is the only nonzero element. If B tj is such an 
element, then B kj = 0 for all k ^ and B i+k j +m = 0 for all 
k > 0 and m< k. This equation relates the 'dependent' vari- 
able S sj to 'independent' variables according to 

(3+df-3d „ 

K = ~ E it**- (15) 



k=j+l 



The number of independent equations is smaller than or equal 
to [(3 + df — 3d]n If the number of independent equations 
is larger than [(3 + d) — 3d], a solution does not exist for S s , 
and the two superspace groups are shown to be inequivalent. 
Alternatively, the number of independent equations can be 
equal to [(3 + df — 3d], then defining a unique solution for S s . 
Finally, the number of independent equations can be smaller 
than [(3 + df — 3d], resulting in more than one solution to 
equation (15). Once values for the independent variables have 
been chosen, equation (15) can be used to compute the 
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remaining variables of S s and S s . For each solution S s of 
equation (15), equation (12) may or may not provide a solu- 
tion for the translational parts of the transformation. 

The strategy for finding the transformation S is now as 
follows. For each trial set of integers for the independent 
variables S sk , check that all dependent variables S • compute to 
have integer values and that det(i?) = 1 and det(e) = ±1. If 
not, discard the trial set. If so, use the values S sj (both 
dependent and independent) in equation (12) and explore 
trial integer sets of variables S vi in search of a modulo 1 
solution. If a solution is found, then equations (11) and (12) 
are both satisfied and the two sets of superspace-group 
operators, G 1 and G 2 , represent distinct but equivalent 
settings of the same superspace group. If no solution is found, 
then we can assume that the groups are not equivalent, 
provided that we have a robust algorithm that searches a 
sufficiently wide range of trial values for each independent 
variable as to guarantee a solution to equations (11) and (12) 
provided one exists. The method of choosing these variable 
exploration ranges is described in the supplementary mate- 
rial. 2 

The number of variables, and therefore the computational 
complexity of the search, increases with the dimension d of the 
modulation. Furthermore, the goal of the proposed analysis is 
to determine which superspace group in the SSG(3 + d)D 
tables a user-given superspace is equivalent to. Since the 
number of superspace groups strongly increases with d, the 
number of candidate equivalencies that need to be tested 
increases dramatically with increasing dimension of the 
superspace, easily reaching several hundreds of groups in the 
worst case (orthorhombic symmetry). Thus we need an algo- 
rithm for evaluating the possible equivalence of two super- 
space groups that is not only robust but also efficient. The 
efficiency of the algorithm boils down to finding the most 
restrictive number of trial sets of integers for which equiva- 
lence needs to be tested (see the supplementary material). 

An algorithm based on these rules has been implemented in 
the software SSG(3 + d)D (Stokes et al, 2011a). For any user- 
given set of superspace operators, the web tool finds sg 
determines the complete list of operators (modulo lattice 
translations) of the superspace group that they generate, as 
well as a minimal list of generators, identifies the equivalent 
superspace group in the SSG(3 + d)D tables and provides the 
coordinate transformation S to the standard setting [equation 
(7)]. 

3. Alternate settings of (3 + 1)-dimensional superspace 
groups 

3.1 . The basic structure space group 

An important reason for the occurrence of non-standard 
settings of superspace groups is the common use of different 
standard settings for superspace groups and three-dimensional 

2 Supplementary material for this paper is available from the IUCr electronic 
archives (Reference: PC5018). Services for accessing this material are 
described at the back of the journal. 



space groups. Structural analysis of modulated crystals often 
proceeds by the initial determination from the main reflec- 
tions of the periodic basic structure along with its three- 
dimensional space group (the basic structure space group, 
BSG). Subsequently, satellite reflections are considered and 
modulation functions and the superspace group are deter- 
mined. For other substances the incommensurate phase is the 
result of a phase transition, so that the three-dimensional 
space group of the unmodulated structure at ambient condi- 
tions is known independently. This space group, or one of its 
subgroups, is preserved as the BSG of the incommensurate 
phase. 3 

In all these cases the BSG is specified before the symmetry 
of the modulation is considered. It is then a matter of chance 
that the superspace group thus obtained will or will not be in 
its standard setting. These points can be illustrated by space 
group No. 62 with standard setting Prima (Hahn, 2002). In this 
setting, TaSe 0 36 Te 2 is modulated with a = (0, 0, er 3 ) and 
cr 3 = 0.6398 [equation (2); van der Lee et al. (1994)], so that 
the (3 + 1)D superspace-group symbol is Pnma(0, 0, y)000, 
which is the standard setting for superspace group No. 62.1.9.1 
in SSG(3 + d)D. 

Thiourea has a lattice-periodic structure with space group 
Prima at ambient conditions. Below T l = 202 K it develops an 
incommensurate modulation with a — (0, a 2 , 0) in the Prima 
setting (Gao et al., 1988; Zuniga et al, 1989). Combining the 
BSG and modulation wavevector leads to the (3 + 1)D 
superspace group Pnma(0, fi, 0)000. SSG(3 + d)D shows that 
this is an alternate setting of superspace group No. 62.1.9.3, 
for which the standard symbol is Pbnm(0, 0, y)000. The 
augmented matrix S that transforms coordinates from the 
original (unprimed) to the standard (primed) settings [equa- 
tion (6)] is given in SSG(3 + d)D as 



/o 


0 


1 


0 


0\ 




(0 


1 


0 


0 


0\ 


1 


0 


0 


0 


0 




0 


0 


1 


0 


0 


0 


1 


0 


0 


0 




1 


0 


0 


0 


0 


0 


0 


0 


1 


0 




0 


0 


0 


1 


0 




0 


0 


0 






\o 


0 


0 


0 


1/ 



(16) 

According to equation (8) the new basis vectors of the basic 
structure are obtained as the upper-left 3x3 part of the 
transpose of S 1-1 . Inspection of equation (16) shows that the 
basis vectors of the basic structure in the standard (primed) 
setting Pbnm(Q, 0, y)000 are obtained by a transformation 
of the basis vectors in the original (unprimed) setting 
Pnma(0, fi, 0)000 as 

a l = a 3> a 2 = a l> a 3 = a 2- (17) 

The fourth row of S shows that the modulation wavevector 
remains the same, but its components with respect to the 
transformed basic structure reciprocal basis vectors are 
obtained by equation (8), 

3 The modulated phase often appears as a twinned crystal if the BSG is a 
subgroup of the three-dimensional space group of the periodic phase (van 
Smaalen, 2007). 
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a' = [(0, 0, 0) + (0, tr 2 , 0)] 



/o o A" 1 

1 0 0 
\0 1 0/ 



2a 3 * + q' = -a x * + q 



(22) 



= (0, 0, cr 2 ) = (0, 0, a 3 ), 



(18) 



in accordance with the standard setting of superspace group 
No. 62.1.9.3. 

3.2. Choice of the modulation wavevector 

A second source of variation of settings is the freedom in 
the choice of the modulation wavevector. Given a modulation 
with modulation wavevector q, any reciprocal vector 



q' = Mja* + « 2 a 2 + « 3 a* ± q, 



(19) 



where n i are integers is an appropriate choice for the modu- 
lation wavevector. A common choice is to select q within the 
first Brillouin zone of the basic structure, i.e. to choose the 
shortest possible vector [equation (19)]. This choice does not 
necessarily correspond to the standard setting of the super- 
space group. 

Transformations that change settings have been extensively 
discussed in van Smaalen (2007). The principles are illustrated 
by the symmetries of A 2 BX 4 ferroelectric compounds with 
the /3-K 2 S0 4 structure type and orthorhombic symmetry 
according to space group No. 62 (Hahn, 2002). Basic structures 
have been described in the standard setting Pnma for some 
compounds, but the most frequently employed settings are 
Pnam and Pmcn (Cummins, 1990). 

K 2 Se0 4 develops an incommensurate modulation below 
T i = 129.5 K with a — (o 1 , 0, 0) in the Pnam setting (Yamada 
& Ikeda, 1984). The incommensurate component is a x = | — S 
with S equal to a small positive number that depends on 
temperature. The superspace group is Pnam(a, 0, 0)0ss. 
SSG(3 + d)D shows that this is an alternate setting of super- 
space group No. 62.1.9.6 Pmcn(Q, 0, y)s00, involving a trans- 
formation of basis vectors and the selection of an alternative 
modulation wavevector according to 



a 3 , a 2 — a 2 , 



-a 1; q=a* + q. 



(20) 



With respect to the transformed reciprocal basis vectors, the 
components of the modulation wavevector are [equation (8)] 



ff' = [(l,0,0) + (a 1 ,0,0)] 



= (o,o,-± + a). 



/0 0 -l\ 
0 1 0 

vi o o ; 



(21) 



The transformed modulation wavevector has a negative 
component and a length larger than a|, which might be 
considered an unfavorable situation. SSG(3 + d)D contains the 
tool transf ormssg, with which any user-specified transfor- 
mation can be applied to the reciprocal basis vectors and 
modulation wavevectors. Employing this tool with [equation 
(20)] 



shows that a modulation wavevector with components 
[equation (8)] 



/' = (o,oj + 5 ) 



(23) 



again represents the standard setting of superspace group No. 
62.1.9.6 Pmcn(0, 0, y)s00. Alternatively, the transformation 



q 

a" 



= -«3 - q = -q 

=K-«) 



(24) 



leads to a non-standard setting of superspace group No. 
62.1.9.6. 

The analysis of symmetry alone does not consider numerical 
values of lattice parameters or modulation wavevectors. 
Therefore, findssg does not employ this information. 
Accordingly, it is impossible to give preference to one of the 
transformations of equation (20) or equation (22). Instead, the 
tool transf ormssg can be used for transformation to the 
desired values. 

Rb 2 ZnCl 4 is incommensurately modulated below T i = 
375 K with a — (0, 0, er 3 ) in the Pmcn setting (Hogervorst, 
1986). The incommensurate component is er 3 = \— S with 
S = 0.04 at room temperature. The superspace group is 
Pmcn(0, 0, y)ss0. SSG(3 + d)D shows that this is another 
alternate setting of superspace group No. 62.1.9.6 
Pmcn(0, 0, y)s00. The transformation now only involves the 
choice of a different modulation wavevector: 



3j — a 1; a 2 — a 2 , 
a = (0, 0, 1 + ct 3 ) = { 0, 0 



a 3 = a 3 , q=a 3 * + q 
4 
3 



(25) 



Like in the previous example, the transformation given by 
findssg does not lead to the setting with the shortest possible 
modulation wavevector for the case of Rb 2 ZnCl 4 . Employing 
transf ormssg shows that the standard setting of superspace 
group No. 62.1.9.6 can also be obtained by the transformation 
of modulation wavevector 



q' = 2a % - q = a£ - q 

a" = (0, 0, 1 - a 3 ) = (o,oJ + s). 



(26) 



As discussed in van Smaalen (2007), replacement of q = er 3 af 
by q = (1 — cr 3 )a| may change the apparently intrinsic trans- 
lational component along the fourth coordinate for symmetry 
operators that possess a nonzero intrinsic translational 
component in the direction corresponding to the incommen- 
surate component of the modulation wavevector. In the 
present example that is (c, s) [mirror operation with intrinsic 
translation (0, 0, \ , |)] being replaced by (c, 0) [mirror opera- 
tion with intrinsic translation (0, 0, 1 , 0)]. These two settings of 
the superspace group correlate with different normal-mode 
descriptions of the same phase transition, for which it has been 
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Table 1 

Symmetry operators of superspace group No. 12.1.8.5 in mixed, BSG and supercentered 
settings. 

The second column refers to the mixed setting chosen by Schutte & de Boer (1993) [who 
incorrectly give the symbol C for the centering translation (| , \ , 0, j)] with the corresponding 
supercentered setting in the third column, featuring the tentative symbols C and C' c for the 
centerings according to Table 3.9 in van Smaalen (2007). The fourth and fifth columns give the 
standard BSG and standard supercentered settings as provided by SSG(3 + d)D. The notation 
of symmetry operators follows SSG(3 + d)D, where — x has been replaced by x and similar. 



Symmetry 
operator 


C'2/m(0P !>0 


C' c 2/m(0p0)s0 


B2/m(0 1 )/)00 


B b 2/m(Q0y)Q0 


Centering 


f 1 1 




(I. 5-0,1) 


(l.o, 


|,o) 




(1,0,1,0) 








(0,0,1,1) 








(0,1,0,1) 








(b bbv 








(1111) 

V2 ' 2 > 2 > V 


Identity 


(x,y, 


z,i) 


(X, Y, Z, T) 


(x,y. 


z,t) 




(X, Y, Z, T) 


Twofold rotation 


(x,y, 


z, z + 1 + D 


(X, Y,Z,T + 1) 


(*, y, 


z,y + 


i) 


(x, Y, z, r) 


Inversion 


(x, y, 


z,i) 


(X, 7, Z, 7) 


(*. y, 


z,t) 




(X, 7, Z, 7) 


Mirror 


(*. y, 


z, z + r + D 


(X, Y,Z,T + 1) 


(x,y, 


z,y + 


t) 


(X, Y, Z, 7) 



setting as opposed to the BSG setting, where a 
modulation wavevector with rational components 
is combined with the standard centered setting of 
the BSG (Stokes et al, 2011a). SSG(3 + d)D 
provides symmetry operators for both the BSG 
and supercentered settings (Stokes et al, 20116). 

As an example, consider blue bronze 
K 03 MoO 3 . (Despite the decimal subscript in the 
usual form of the chemical formula, it possesses a 
fully ordered crystal structure with two formula 
units K 3 Mo 10 O 30 per unit cell.) Blue bronze 
develops an incommensurate charge-density 



wave (CDW) below T c 
modulation wavevector 



183 K with a 



describes a distortion in terms of a A 2 soft optical phonon, 
while 03 = 5 + S leads to the preferred description of the 
distortion in terms of a A 3 soft optical phonon (Axe et al, 
1986). It is well known that a change of setting will sometimes 
change the irreducible representation that contributes to a 
distortion without changing the physical distortion itself. 

3.3. The supercentered setting 

Aperiodic crystals are characterized by d modulation 
wavevectors, each of which possesses at least one irrational 
component. According to symmetry, the values of these 
components can be viewed as variable rather than having 
specific irrational values. The other components are either 
zero or may assume rational values as allowed by the point 
symmetry. It is easily checked against the list of Bravais classes 
of (3 + d)D superspace groups that the allowed rational 
components of modulation wavevectors are \ and \ in the 
standard superspace-group settings, as well as nonzero inte- 
gers in the case of BSGs based on a centered lattice of the 
basic structure. The modulation wavevector is usually sepa- 
rated into a rational part, with zeros and rational numbers as 
components, and an irrational part, with zeros and the variable 
components, according to 



q — q r + 4. a = a r + a 1 . 



(27) 



A modulation wavevector with nonzero rational compo- 
nents may naturally occur when a diffraction pattern of a 
modulated crystal is indexed, first determining the unit cell - 
centered if required - of the basic lattice, and then selecting an 
appropriate modulation wavevector, e.g. the shortest possible 
vector. However, rational components of q imply that the 
lattice in superspace is a centered lattice with the special 
feature that centering vectors contain nonzero components 
both along the fourth coordinate and along at least one of the 
three physical coordinates. Employing centered unit cells for 
centered lattices is common practice. It has several advantages 
in crystallographic analysis, facilitating the description of 
reflection conditions and the analysis of point symmetry. For 
superspace groups it has been denoted as the supercentered 



a = (0, 0.748, 0 



(28) 



at T = 100 K in the setting of the high-temperature space 
group. Schutte & de Boer (1993) have determined the crystal 
structure of the incommensurate phase. With the modulation 
wavevector of equation (28), they obtained the (3 + 1)D 
superspace group C'2/m(0, j3, £)s0. This mixed setting contains 
the C center (| , \ , 0, |), which possesses a nonzero component 
along the fourth coordinate (van Smaalen, 2007), and which is 
different from the BSG setting comprising the C center 
(I , \ , 0, 0). f indssg shows that C'2/m(0, fi, ^)s0 is an alternate 
setting of superspace group No. 12.1.8.5 B2/m(0, \ , y)00 
(Table 1 and Fig. 1). In physical space, the transformation from 
the setting of Schutte & de Boer (1993) to the standard BSG 
setting involves a permutation of the unit-cell axes and the 
choice of a new modulation wavevector according to 



a 7 = — i 



3-> 



-q 



(29) 



Furthermore, an origin shift of (| , \ , 0, 0) is required in order 
to bring the origin onto the operator (2, 0) instead of (2, s). 
The components of q with respect to the transformed (primed) 
reciprocal basis vectors follow from equation (8) as [compare 
equation (21)] 



(0,1,0)- 0,<r 2 



1 



(\ 0 0\ 

0 0 1 
\0 -1 0/ 



(0,i,l-a 2 ) = (0,i,a 3 ). 



(30) 



The standard supercentered setting is obtained from the 
standard BSG setting B2/m(0\a 3 )0Q by the transformation of 
superspace basis vectors as given in SSG(3 + d)D, and corre- 
sponds to the following transformation of the physical-space 
basis vectors and modulation wavevector (capital letters 
indicate the supercentered setting): 



A l= a l> A 2 = 

ct 1 = (0, 0, a 3 ). 



2a 2 , 



(31) 



Observe that the modulation wavevector is purely irrational in 
the supercentered setting, as expected [equation (27)]. 
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findssg 

Input setting 

Centering 

(0,0,0,0); (1/2,1/2,0,1/2) 
Operators 

(-x,y,-z,-z+t+l/2); (-x,-y,-z,-t); (x,y,z,t); (x,-y,z,z-t+l/2) 

Standard settings 

Superspaee group: 12.1.8.5 B2/m(0,l/2,g)00 [Y:1.52] 
Bravais class: 1.8 B2/m(0,l/2,g) [JJdW:1.8] 

Transformation to supercentered setting: Al=al, A2=2a2+a4, A3=a3, A4=a4 

BASIC SPACE GROUP SETTING 

Modulation vectors: ql'=(0,l/2,g) 

Centering: (0,0,0,0); (1/2,0,1/2,0) 

Non-lattice generators: (-x,-y,z,-y+t); (x,y,-z,y-t) 

Non-lattice operators: (x,y,z,t); (-x,-y,z,-y+t); (-x,-y,-z,-t); (x,y,-z,y-t) 

SUPERCENTERED SETTING 

Modulation vectors: Ql -(0,0,G), where G=g 

Centering: (0,0,0,0); (1/2,0,1/2,0); (0,1/2,0,1/2); (1/2,1/2,1/2,1/2) 

Non-lattice generators: (-X,-Y,Z,T); (X,Y,-Z,-T) 

Non-lattice operators: (X,Y,Z,T); ( X,-Y,Z,T); (-X,-Y,-Z,-T); (X,Y,-Z,-T) 
Reflection conditions: HKLM:H+L=2n; HKLM:K+M=2n 

Affine transformation to standard basic space group setting 

S * g(input) * S - g(standard) 

where g is an augmented matrix for an operation in the superspaee group. 

Also, S * r(input) - r(standard), 

where r is an augmented position vector, (x,y,z,t,l). 

10 0 0 1/2 10 0 0 -1/2 

0 0 -10 1/2 0 0 10 0 

S = 0 1 0 0 0 S"' = 0 -1 0 0 1/2 

0 10 -10 0 01-1 0 

00001 0000 1 



al'-al 
a2' = -a3 
a3' = a2 

al =al' 
a2 = a3' 
a3 = -a2' 



groups, where the latter refer to 
incommensurate crystals with two 
independent modulation waves, while 
the (3 + 2)D superspaee groups refer 
to crystals with two symmetry- 
related modulation waves, such as 
those in Bravais class No. 2.57 
P4/m(a, 0, 0)(-p, a, 0). 

As an example of a (3 + 1 + 1)D 
superspaee group, consider 

P2/m(0, 0, Xi>0(0, 0, y 2 )s0. findssg 
shows that this is an alternate setting 
of superspaee group No. 10.2.5.6 
P2/m(0, 0, Xi>0(0, 0, y 2 )00. The trans- 
formation which brings the original 
setting into the standard BSG setting is 



(33) 



al»' = al» 
a2*' = -a3* 
a3»' = a2» 

al» = al*' 
a2* = a3*' 
a3» = -a2«" 

q l ' = -q l +a2* = (0,l/2,g) 
ql =-qr + a3*' = (0,-g+l,l/2) 



Figure 1 

Output of findssg on SSG(3 + d)D, showing the equivalence of superspaee group C'2/m(0o 2 |)s0 to 
superspaee group No. 12.1.8.5 52/m(0 |cr 3 )00. 



4. A plethora of settings 

4.1. (3 + 2)-Dimensional superspaee groups 

4.1 .1 . General features. Different settings of (3 + d)D 
superspaee groups are obtained when different settings of the 
BSG are chosen (§3.1). For d = 1, further settings result from 
the ambiguity in the choice of the modulation wavevector: all 
modulation wavevectors that differ by a reciprocal-lattice 
vector of the basic structure are equally valid [equation (19)] 
and may define different settings of a superspaee group (§3.2). 
For d>2, additional coordinate transformations involve the 
replacement of the modulation wavevectors by linear combi- 
nations of them. 

Generalizing equation (7) shows that any set of reciprocal 
vectors (;' = 1, . . . , d) 



The difference between the two settings 
is that the operator (2, ss) in one setting 
is equivalent to the operator (2, sO) in 
the standard setting. The transforma- 
tion modifies the intrinsic translation of 
a twofold axis along one of the super- 
space coordinates. This is a feature 
specific to transformations of the type 
of equation (33), while the modification 
of the modulation wavevector by a 
reciprocal-lattice vector of the basic 
structure [equation (19)] can only affect 
the intrinsic translations of symmetry 
operators that are screw axes or glide planes in three dimen- 
sions. More complicated linear combinations of modulation 
wavevectors may be required, as in the transformation 
between 74(0, 0, y t )q(0, 0, y 2 )s and the standard setting 
74(0, 0, ftMO, 0, y 2 )0 of superspaee group No. 79.2.62.3 
[equation (32)]: 



2q X 



(34) 



(32) 



/=i 



The same concept can be applied for reducing the 
number of rational components of the modulation 
wavevectors. SSG(3 + d)D shows that superspaee group 
P222(|, p u 0)00(|,^ 2 , 0)00 is an alternate setting of super- 
space group No. 16.2.19.3 P222(0, \ , K,)00(0, 0, y 2 )00 with 
y 1 = p i and y 2 = ~P l + P 2 . The transformation between 
these settings involves a linear combination of the two 
modulation wavevectors as well as a change of the setting of 
the BSG according to 



is an appropriate choice for the set of d modulation wave- 
vectors, where S Mji and S S j k are integers and det(5 , E/Jt ) = ±1. 
Such linear combinations of the modulation vectors of the 
standard setting are often necessary to make the description of 
experimental diffraction data simpler and more intuitive, but 
can also have the opposite effect if applied arbitrarily. 

In analyzing superspaee groups with d = 2, it is useful to 
distinguish between (3 + 2)D and (3 + 1 + 1)D superspaee 



q =q\ 



a 2 — a l > 

,2 

q = - 



(35) 



where again the primed vectors refer to the standard BSG 
setting. 

4.1.2. NbSe 3 . Several of the features discussed here are 
illustrated by the example of NbSe 3 . NbSe 3 develops an 
incommensurate CDW below r CDW1 = 145 K. A second, 
independent CDW develops below r CDW2 = 59 K, then 
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resulting in an incommensurately 
modulated structure with two inde- 
pendent modulation waves, and 
with symmetry given by the 
(3 + 2)D superspace group No. 
11.2.6.4 P2Jm{\ , 0, ft)00(0, 0, y 2 )00 
from SSG(3 + d)D. Inspection of the 
list of superspace groups shows that 
No. 11.2.6.4 is the only superspace 
group in its Bravais class that has BSG 
P2 1 /m. This implies that any possible 
combinations of nonzero intrinsic 
translations along the fourth and fifth 
coordinate axes are equivalent to the 
setting FIJmQ , 0, fr)00(0, 0, y 2 )00 by 
a suitable transformation. 

The modulated, low-temperature 
crystal structure of NbSe 3 has been described in superspace 
group P2 1 /m(0, 0)s0(|, /J 2 , |)00 (unique a 2 axis) with (van 
Smaalen et al, 1992) 

'0 




Figure 2 

Projection of the basic structure of NbSe 3 . All atoms are in mirror planes. Hatched and open circles 
are atoms at 0.25 and 0.75 of the projected coordinate, respectively. Small circles are Se; large circles 
are Nb, with numbers indicating Nbl, Nb2 and Nb3 atoms. Symmetry operators (2j, 00) and (2 l , sO) 
alternate in the supercentered setting. Unit cells are indicated for the published BSG setting (a x , a 3 



solid lines), the standard BSG setting [a\ 



supercentered setting [A ls A 2 = a' 2 ; dashed lines; see equation (41)]. 



a' 2 ; dashed lines; see equation (37)] and the 



_ /0 cr 12 0\ _ /0 0.241 0\ 

U ^ \)~\\ °- 260 V 



(36) 



This setting naturally arises for the following choices: 

(i) The BSG is equal to the space group of the periodic 
structure at ambient conditions, which has a unique axis a 2 
that is the preferred setting for monoclinic three-dimensional 
space groups. 

(ii) The choice of axes a t and a 3 is that of the previously 
determined periodic crystal structure at ambient conditions. 

(iii) Modulation wavevectors are chosen within the first 
Brillouin zone. 

(iv) The first modulation wavevector is that of the first 
CDW and the second modulation wavevector applies to the 
second CDW. 

All four choices need to be adapted, in order to arrive at the 
standard setting of this superspace group: 

(i) The standard setting of the superspace group has 
incommensurate components of the modulation wavevectors 
along a*, thus requiring a reordering of the basic structure 
axes. 

(ii) Transforming the second modulation wavevector into a 
wavevector with one nonzero rational component requires a 
basic structure monoclinic unit cell that involves linear 
combinations of the axes aj and a 3 [compare to equation (1)]. 
Notice that this transformation does not affect the symbol of 
the BSG. 

(iii) The transformation of (2 1 , 5-0) into (2 1 , 00) requires the 
transformation q' 1 = a| ± q 1 [equation (32)]. 

(iv) The standard setting requires interchanging the two 
modulation wavevectors. 

Altogether, the transformation from the published setting 
to the standard BSG setting of superspace group No. 11.2.6.4 
is achieved by 



a i — a^> 



= a, 



q = q^ 



= a 2 * + q\ 



(37) 



which implies a transformation of reciprocal basis vectors of 
the basic structure as 



a'* 



= a, + a 3 



i* 
a 2 



(38) 



The components of the modulation wavevectors with respect 
to the transformed reciprocal basis vectors follow from 
equation (8) or by inspection of equations (37) and (38): 



,(\ 0 0 a 22 \ 

\0 0 o> 2 J \Q 0 l+a 12y l 



(39) 



The tool transforms sg can be used to demonstrate that an 
alternate transformation, defined by a different choice of the 
second modulation wavevector, also leads to the standard 
BSG setting of superspace group No. 11.2.6.4: 



a": 


= a' 1 = 


a 3 , 


a" 2 


= a' 2 


— a x a 3 










= q 4 




a" 3 


= a 3 = 


a 2 , 


q" 1 


= q 2 


q" 2 


= 2a'* - 


a 

- q 


= a*- 


q 1 





\0 0 &') Vo 0 \-oJ 



(40) 



Choices (i), (ii) and (iv) are arbitrary - there does not appear 
to be a compelling reason to adhere to the standard setting 
except to establish the equivalence of different crystal struc- 
tures. The choice (iii) of the modulation wavevector is related 
to the important question about the real wavevectors of the 
CDWs, which is not obvious because the incommensurate 
components of the modulation wavevectors can either be 
~0.25 or ~0.75, depending on the setting. This is most easily 
analyzed with the help of the supercentered setting, which 
follows from the standard BSG setting by the transformation 



2a"! , A 2 
0 
0 



/0 0 a" 13 \ 
V0 0 a" 2 J" 



(41) 



Structural analysis has shown that the first CDW (q" in the 
standard setting) is located on a pair of chains of niobium 
atoms, denoted as the Nb3 atoms, while the second CDW (q" 1 
= q' 1 ) is located on a pair of chains of Nbl atoms (Fig. 2). 
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Table 2 

Superspace groups for incommensurate compounds with two-dimensional modulations. 



Given are the published modulation wavevectors and superspace-group symbols, the number and symbol of the standard BSG setting of the superspace group in 
SSG(3 + d)D, and the transformation of the published basic structure unit cell to the standard BSG setting as well as the modulation wavevectors of the latter. 









Published 


SSG(3 + d)D 


Superspace-group symbols for 


Standard BSG setting 


Compound 


Note 


r(K) 


q'/q 2 


No. 


published/standard BSG settings 


a', 


a 2 ,a 3 


q'Vq' 2 


Mo 2 S 3 


(«) 


295 


(-0.056, 0.5, 0.229) 


2.2.1.1 


Pl{a, p, y)(k, ix, v)T 


a l 


a 2 , a 3 


(-0.056, 0.5, 0.229) 








(0.5,0.469,0) 




-Pl(<*i, A, fr)0(<x 2 , ft, K 2 )0 






(0.5, 0.469, 0) 


(Bi,Pb) 2 (Sr,Bi,Pb,Ca) 2 - 


(b) 


295 


(0.234, 0, 0) 


9.2.4.1 


P -.A'a : PI 


a 3 


a,,a 2 


(1, 0.234, 0) 


CuO„ + , 






(0.144, 0, 0) 




B6(o! 1 ,A,0)0(a 2>/ 8 2 ,0)0 






(0, 0.144, 0) 


NbSe 3 


(c) 


<59.0 


(0, 0.241, 0) 


11.2.6.4 


P : P2 t /m(0, p, 0)(1 , 0, 1>1; 11 


% 


a l + a 3. a 2 


(1 , 0, 0.740) 








(1,0.260,1) 




P2 1 /m(\ , 0, tt)00(0, 0, y 2 )00 






(0, 0, 0.759) 


TTF TCNQ 


(d) 


<38.0 


(0.25, 0.295, 0) 


14.2.16.6 


P : P2 l /c : cmm 


"3 


a l > a 2 


(0, 0.25, 0.295) 








(-0.25, 0.295, 0) 




P2 1 /c(a, p, K)00(-ff, -p, y)00 






(0, -0.25, 0.295) 


(P0 2 ) 4 (W0 3 ) 8 


(e) 


<80.0 


(0.330, 0.292, 0) 


19.2.50.3 


P2 1 2 1 2 1 (a, p, 0)(-a, p, 0) 


a 3 


aj, a 2 


(0, 0.330, 0.292) 








(-0.330, 0.292, 0) 




P2 1 2 1 2 1 (0, p, y)000(0, -p, y)000 






(0, -0.330, 0.292) 


S m.2/3 S4 


if) 


295 


(0.3333, 0.3217, 0) 


62.2.50.22 


Pmnb(a p 0)gmg 


a 3 


a 2 , -a! 


(0, 0.6783, 0.3333) 






(-0.3333, 0.3217, 0) 




Pbnm(Q, p, k)000(0, ~p, y)000 






(0, -0.6783, 0.3333) 


GdS lg2 


(g) 


293.0 


(0.25,0.33,1) 


85.2.58.2 


PA/n{a, p, j)(00)(ii) 


a i 


a 2 , a 3 


(0.25,0.33,1) 








(-0.33,0.25,1) 




PA/n(a, p, 1)00(-^, a, 1)00 






(-0.33,0.25,1) 


LaSe 185 


(h) 


130.0 


(-0.292, 0.292, 1) 


85.2.58.2 


P4/n(-a, p, l)(/3, a, 1)00 


a i 


a 2 , a 3 


(0.292, 0.292, 1) 








(0.292, 0.292, 1) 




P4/n(a, p,l)00(-p,a,l)00 






(-0.292, 0.292, 1) 


Bao.3 9 Sr 0 . 61 Nb 2 0 6 


(0 


295 


(0.3075, 0.3075, 1) 


100.2.69.13 


P4bm(p pi, p -p 1)000 


»i 


a 2 , a 3 


(0.3075, 0.3075, 1) 








(0.3075, -0.3075, 1) 




P4bm(a, a, l)000(-a, a, 1)000 






(-0.3075,0.3075,1) 


Sr 2 TiSi 2 0 8 


0') 


295 


(0.3, 0.3, 1) 


100.2.69.14 


P4bm(—a, a,\; a, a, |)0gg 


»i 


a 2 , a 3 


(0.3, 0.3, 1) 








(-0.3,0.3,1) 




P4bm(a, a, l)00s(-a, a, 1)000 






(-0.3,0.3,1) 


CajCoSijO, 


(*) 


295 


(0.2913, 0.2913, 0) 


113.2.68.6 


P42 x m{p p 0, — p p 0)4mg 


a i 


a 2 , a 3 


(0.2913, 0.2913, 0) 








(-0.2913, 0.2913, 0) 




P42 1 m(a, a, 0)00s(-a, a, 0)000 






(-0.2913, 0.2913, 0) 


(Sr 0 13 Ca 087 ) 2 CoSi 2 O 7 


(0 


295 


(0.286, 0.286, 0) 


113.2.68.6 


P42 x m: 4mg 


a i 


a 2 , a 3 


(0.286, 0.286, 0) 








(-0.286, 0.286, 0) 




P42 1 m(a, a, 0)00s(-a, a, 0)000 






(-0.286, 0.286, 0) 


CaNdGa 3 0 7 


(m) 


293 


(0.2319, 0.2319, 0) 


113.2.68.6 


P42 x m(a, a, 0)00s(-a, a, 0)000 


a i 


a 2 , a 3 


(0.2319, 0.2319, 0) 








(-0.2319, 0.2319, 0) 




P42 1 m(a, a, 0)00s(-a, a, 0)000 






(-0.2319, 0.2319, 0) 


CaLaGa 3 0 7 


(n) 


293 


(0.2160, 0.2160, 0) 


113.2.68.6 


P42 lt n(a, a, 0)00s(-a, a, 0)000 


a i 


a 2 , a 3 


(0.2160, 0.2160, 0) 








(-0.2160, 0.2160, 0) 




P42 1 m(a, a, 0)00s(-a, a, 0)000 






(— O.ztoU, 0.2160, 0) 


Ni 581 SnTe 2 


(0) 


153 


(0, -0.410, 0) 


139.2.67.7 


I4/mmm(0 — a 0, a 0 0)0. ss. mm 


a i, a 3 


(0.410, 0, 0) 








(0.410, 0, 0) 




14/mmm(a, 0, 0)0s00(0, a, 0)0s00 






(0, 0.410, 0) 


ir-TaS 2 


ip) 


< 350.0 


(0.245, 0.068, 1) 


147.2.72.1 


P3(a p 1/3) 


a,, a 2 , a 3 


(0.068, -0.313,1) 








(-0.068, 0.313, D 




P3(a, ,S, l)00(-a - /3, a?, |)00 






(0.245, 0.068, 1) 


l'-Cu 3+I Si 


(?) 


298 


(0.244, 0.244, 1) 


162.2.76.3 


raim(a, a, l)00(-2a, a, 1)00 


a i 


a 2. % 


(0.244, 0.244, 1) 








(-0.488, 0.244, 1) 




P31m(a, a, l)00(-2a, a, 1)00 






(-0.488, 0.244, 1) 


4H b -TaSe 2 


w 


< 410.0 




176.2.80.1 


F6 3 /m(a 0 0) 


a i 


a 2 , a 3 


(2- 1- 01 

M3 ' 13 ' ™> 








(-£.£.0) 




F6 3 /m(a, /3, 0)00(-a - /S, a, 0)00 






(- A ^ 0) 

V 13 ' 13 ' u ' 


Ag,,. 6 NbS 2 


(*) 


100 


(0.175, 0.175, 0) 


186.2.83.4 


P6 3 mc : 6mm 


a 2 , a 1? a 3 


(0.175,0.175,0) 








(-0.175, 0.350, 0) 




P6 s mc(a, a, 0)000(-2a, a, 0)000 






(-0.350, 0.175, 0) 




(0 


295 


(0.43,0.43,0) 


194.2.83.4 


P6^/mmc : 6\mm 


a 2 


»1. ~ a 3 


(0.43, 0.43, 0) 








(-0.43,0.86,0) 




P6 3 /mmc(a, a, 0)0000(-2a, a, 0)0000 






(-0.86, 0.43, 0) 


Cu 5 Sb 


('<) 


295 


(0.43,0.43,0) 


194.2.83.4 


P6^/mmc : 6lmm 


a 2 


■l. ~ a 3 


(0.43, 0.43, 0) 








(-0.43,0.86,0) 




P6 3 /mmc(a, a, 0)0000(-2a, a, 0)0000 






(-0.86, 0.43, 0) 



Notes: (a) Schutte et al. (1993). (b) Gao a al (1990). (c) van Smaalen a al (1992). (d) TTF (tetrathiafulvalene) and TCNQ (tetracyanoquinodimethane), Bouveret & Megtert (1989). (e) 
Ludecke el al (2001); an extensive review of basic structures and CDW transitions of the phosphate bronzes (P0 2 ) 4 (W0 3 ) 2m (2 < m < 14) is given by Roussel et al (2000). (f) Lafond et 
al. (1996). (g) Tamazyan et al. (2003). (h) Graf & Doert (2009). (i) Woike et al. (2003). (j) Hoche et al. (2003). (k) Hagiya et al (1993). (0 Bagautdinov et al. (2000). (m) Wei et al. (2011). 
(«) Wei et al. (2012). (o) Isaeva et al (2007). (p) Yamamoto et al (1990) and Spijkerman et al. (1997). (q) Palatinus et al. (2011). (r) Ludecke et al. (1999). (s) van der Lee et al (1991). (i) 
Yamamoto (1983). (u) Motai et al. (1993). 



SSG(3 + d)D provides the explicit form of the symmetry 
operators in the supercentered setting. Employing these 
operators, one finds that the double chain of Nbl atoms 
centered on , \ , Z, T, U) of the supercentered unit cell is 
located on the screw axis [2 Z 11\ \ , 1, \ , \ , 0}. 4 This is a screw 
operator (2 X , sO) as it is generated by the combination of the 
screw (2 U 00) and the centering translation {J , 1, 0, \ , 0}. The 
pair of chains of Nb3 atoms is related by the operator 
{2 Z 11|0, 1, \, 0, 0}, which is a screw operator (2 l , 00). We judge 



4 The symbol {2 2 11 1 1 , 1, \ , \ , 0} indicates an operator with a twofold rotation 
2 Z of the BSG along a 3 combined with a unit 2x2 matrix s, and followed by a 
superspace translation (1,1,1,1,0). 



that q"? = q" 2 = (0, 0, 0.759) is the real wavevector of the 
Nb3 modulation, because an additional phase shift is not 
involved on application of this symmetry. On the other 
hand, the second wave with q'; = (0, 0, 0.260) [equations (39) 
and (40)] implies symmetry for the pair of Nbl chains 
involving a phase shift of one half. The real wavevector 
thus is af — ql = (0, 0, 0.740), resulting in the setting 
P2 1 /mQ , 0, yi)s0(0, 0, y 2 )00 of superspace group No. 
11.2.6.4. SSG(3 + d)D shows that the standard setting can be 
restored by a shift of the origin of \ along a 1 . With this final 
transformation, the symmetry of NbSe 3 is described in the 
standard setting, and the components of the modulation 
wavevectors show that both CDWs are waves with wave- 
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vectors of ~0.75 on their respective double chains of niobium 
atoms. 

4.1.3. Centerings in internal space. Table 2 compiles 
superspace groups for a series of compounds with two- 
dimensional modulations. Symbols for the superspace groups 
from the original publications encompass a disparate set of 
notations, including symbols based on the online database 
of (3 + d)D superspace groups (d=l,2, 3) of Yamamoto 
(2005), as in the case of Ca 2 CoSi 2 0 7 , symbols based on Janner 
etal. (1983), as in the case of Mo 2 S 3 , and symbols derived from 
these notations, such as replacing p by a in the case of 
Sm 2/3 Cr 2 S 4 , as well as other ad hoc symbols. 

While for several compounds a permutation is required 
of the basis vectors of the basic structure unit cell in 
order to transform the published setting into the 
standard setting, other, less trivial transformations 
occur too. SSG(3 + d)D shows that the symmetry of 
tetrathiafulvalene tetracyanoquinodimethane (TTF TCNQ), 
P2 1 /c(a, fi, y)00(— a, —ft, y)00, is based on a supercentered 
lattice, where the centering exclusively involves the two 
internal coordinates. The supercentered setting has modula- 
tion wavevectors Q 1 = (a u , er 12 , 0) and Q 2 = (0, 0, er 23 ) with 

q 1 = Q 1 + Q 2 = (a n ,a 12 ,a 23 ) 

q 2 = -Q 1 +Q 2 = (-ff 11> -ff 12 , 0 - 23 ), (42) 

resulting in the reflection condition and corresponding 
centering translation: 

(HKLM 1 M 2 ): M l +M 2 = 2n (n is an integer) 
(0,0,0,1,1). ' C43) 

Monoclinic symmetry is the lowest symmetry where 
this kind of superspace lattice centering can occur. For 
Ca 2 CoSi 2 0 7 , on the other hand, the superspace group 
P42 1 m(a, a, 0)00s(— a, a, 0)000 does not possess a centered 
lattice, despite the seemingly simpler modulation wavevectors 
Q 1 = (cr, 0, 0) and Q 2 = (0, a, 0), which are related to q 1 and q 2 
as in equation (42). The reason is that Q 1 and Q 2 are related 
by symmetry in the same way as q 1 and q 2 are, and the would- 
be supercentering does not have an advantage over the 
primitive lattice from the point of view of symmetry. Examples 
of supercentered lattices with higher symmetries are given for 
three-dimensional modulations in §4.2. 

A peculiar feature of the modulation of TTF TCNQ is that 
one of the unrestricted components is not experimentally 
distinguishable from zero (a n = a 21 = aj 3 = a 23 = a — 0). 
The explanation probably lies in the optimal phase relations 
between the CDWs on neighboring stacks of TTF or TCNQ 
molecules, as it is governed by the physics of CDW formation. 
However, in this case the phase relation is not reflected in 
the symmetry of the crystal structure. Similar observations 
can be made for Mo 2 S 3 = a 2 — \ and y 2 = 0), 
(Bi,Pb) 2 (Sr,Bi,Pb,Ca) 2 Cu0 6+5 (a 2 = 0) and LaSe 185 (a = fi). 
For these compounds, the special values of the components 
of the modulation wavevectors are reminiscent of the 
higher symmetries at high temperatures [monoclinic for 
Mo 2 S 3 and orthorhombic for the high-T^ superconductor 



Table 3 

Selected (3 + 2)D superspace groups with acentric trigonal symmetry. 
Other superspace groups exist that differ in the intrinsic translations. 
No. Superspace-group symbol 

143.2.72.1 P3(a, ft, \)0(-a - /3, a, 1)0 

149.2.76.3 P312(a, a, i)000(-2a, a, 1)000 
150.2.78.1 P321(a, 0, l)000(-a, a, 1)000 

143.2.80.4 P3(a, p, 0)0(-a - P, a, 0)0 

149.2.82.6 P312(a, 0, 0)000(-a, a, 0)000 

150.2.82.4 P321(a, 0, 0)000(-a, a, 0)000 

149.2.83.7 P312(a, a, 0)000(-2a, a, 0)000 

150.2.83.5 P321(a, a, 0)000(-2a, a, 0)000 



(Bi,Pb) 2 (Sr,Bi,Pb,Ca) 2 Cu0 6+a ] or the higher symmetry of a 
hypothetical basic structure (Laue symmetry 4/mmm for 
LaSej 85 ). 

The lattice type (primitive, centered BSG or supercentered) 
is the same for all superspace groups belonging to a Bravais 
class. Likewise, the choice of modulation wavevectors should 
be the same for all superspace groups within a single Bravais 
class: the Bravais class is defined by the point symmetry of the 
lattice together with the modulation wavevectors. A further 
requirement on the modulation wavevectors is that they must 
transform according to the three-dimensional point symmetry 
of the superspace group. These requirements become impor- 
tant for the selection of modulation wavevectors in the case of 
trigonal and hexagonal Bravais classes of (3 + 2)D superspace 
groups. Since all these Bravais classes contain superspace 
groups with acentric trigonal symmetry (Table 3), it is neces- 
sary to choose a pair of modulation wavevectors that enclose 
an angle of 120° and not 60° (Fig. 3). With the exception of the 
recent study on r]'-Cu 3+x Si, this condition has not been obeyed 
in studies of the compounds with trigonal or hexagonal 
symmetries listed in Table 2, where the angle between q 1 and 
q 2 was chosen as 60°. While not wrong in these cases, it is 
highly preferable to describe these structures using an angle of 
120° between the modulation wavevectors so as to be 
consistent with the settings of their Bravais classes. 



4.2. (3 + 3)-Dimensional superspace groups 

Different settings of (3 + 3)D superspace groups are 
obtained by means of the same degrees of freedom that apply 
to (3 + 2)D superspace groups. That is, the setting of a 
(3 + 3)D superspace group depends on the choice of basic 
structure basis vectors and on the freedom in the choice of 
modulation wavevectors, including the possibility to replace 
the modulation wavevectors by linear combinations of them 
[equation (32)]. 

4.2.1. Supercentered setting of (TaSe 4 ) 2 l. (TaSe 4 ) 2 I has a 
periodic structure with space group 7422 at ambient condi- 
tions. A CDW develops below T CDW = 263 K. It is expressed in 
the diffraction by the presence of eight incommensurate 
satellite reflections around each main reflection, which can be 
indexed as first-order satellite reflections according to the four 
modulation wavevectors 
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Table 4 

Superspace groups for incommensurate compounds with three-dimensional modulations. 



Given are the published modulation wavevectors and superspace-group symbol, the number and symbol of the standard BSG setting of the superspace group in 
SSG(3 + d)D and the modulation wavevectors in the standard BSG setting. The published and standard basic structure unit cells are equal to each other: a',- = a, 
(i = 1, 2, 3). 











SSG(3 + d)D 


Superspace-group symbols for 




Standard BSG setting 


Compound 


Note 


r(K) 


Published qVqVq 3 


No. 


published/standard BSG settings 




q'Vq'Vq' 3 


(TaSe 4 ) 2 l 


(«) 


< 263.0 


(0.0641,0.0641,0.1510) 


97.3.179.24 


/422(a, a, y)000 




(0.0641,0.0641,0.1510) 








(-0.0641, 0.0641, 0.1510) 




/422(a, a, y)000(-a, a, y)000(a, -a, y)000 


(-0.0641, 0.0641, 0.1510) 








(0.0641, -0.0641,0.1510) 








(0.0641, -0.0641,0.1510) 




(b) 


295 


(0.398, 0, 0) 
(0, 0.398, 0) 
(0, 0, 0.398) 


225.3.209.1 


P : Fm3m : Pm3m 

Fm3m(a, 0, 0)000(0, a, 0)000(0, 0, a)000 




(0.398, 0, 0) 
(0, 0.398, 0) 
(0, 0, 0.398) 


None 


(c) 




(0, P, P) 
08, 0, P) 

08, A o) 


225.3.212.5 


Fm3m(0pp, pOp, pp0)mtm 

Fm3m(0, P, / 6)000( / 6, 0, /S)000(/S, /3, 0)000 




(0,P,P) 

(P,o,p) 

08, P, 0) 


Cu 9 BiS,, 


(d) 


295 


(-0.154,0.154,0.154) 
(0.154, -0.154, 0.154) 
(0.154, 0.154, -0.154) 


225.3.215.7 


P : Fm3m : Fm3m 

f? in \ i j j i cv /v /"i/ImMMi/v /v /v ImOMi /v /v 

1 IllJlllilX, IX, IX JUUUI IX , — IX, — Lt JUUUI — Lt , IX, 




(0.154, 0.154, 0.154) 
(0.154, -0.154, -0.154) 
(-0.154, 0.154, -0.154) 


Bio.85Cr () jsOj 73 


(<0 


295 


(0.2866, 0.2866, 0.2866) 
(0.2866, -0.2866, -0.2866) 
(-0.2866, 0.2866, -0.2866) 


225.3.215.7 


P : Fm3m : Frriim 

Fm3m(a, a, a)000(a, —a, — a)000(— a, a, 


-a)000 


(0.2866, 0.2866, 0.2866) 
(0.2866, -0.2866, -0.2866) 
(-0.2866, 0.2866, -0.2866) 


Bi 084 Mo ()1( ,O lv4 


if) 


295 


(0.2926, 0.2926, 0.2926) 
(0.2926, -0.2926, -0.2926) 
(-0.2926, 0.2926, -0.2926) 


225.3.215.7 


P : Fm3m : Fm3m 

Fm3m(a, a, a)000(a, —a, — a)000(— a, a, 


-a)000 


(0.2926, 0.2926, 0.2926) 
(0.2926, -0.2926, -0.2926) 
(-0.2926, 0.2926, -0.2926) 


Bio.78Nbo.22O1.72 


(g) 


295 


(0.37, 0.37, 0.37) 
(0.37, -0.37, -0.37) 
(-0.37, 0.37, -0.37) 


225.3.215.8 


P : Fm3m : Fd3m 

Fm3m(a, a, a)q00(a, —a, —a)q00(—a, a, 


-a)000 


(0.37, 0.37, 0.37) 
(0.37, -0.37, -0.37) 
(-0.37, 0.37, -0.37) 


^'103^25 0 217 


(h) 


295 


(0.375, 0.375, 0.375) 
(0.375, -0.375, -0.375) 
(-0.375, 0.375, -0.375) 


225.3.215.8 


P : Fm3m : Fd3m 

Fm3m(a, a, a)q00(a, —a, —a)q00(—a, a, 


-a)000 


(0.375, 0.375, 0.375) 
(0.375, -0.375, -0.375) 
(-0.375, 0.375, -0.375) 


BaBi 3 0 55 


(0 


295 


(0, 0.3835, 0.3835) 


229.3.211.5 


P : Im3m : lm3m 




(0, 0.3835, 0.3835) 






(0.3835, 0, 0.3835) 
(0.3835, 0.3835, 0) 




Im3m(0, p, P)000(fi, 0, /3)000(/3, p, 0)000 




(0.3835, 0, 0.3835) 
(0.3835, 0.3835, 0) 


V 6 Ni 16 Si 7 


(7) 


295 


(0.361, 0.361, 1.361) 
(0.361,-0.361,1.361) 
(-0.361, -0.361, 1.361) 


229.3.214.8 


: Im3m : Fm3m 
Im3m(a, a, a)q00(a, —a, — a)q00(— a, a, 


-a)000 


(1.361, 1.361, 1.361) 
(1.361, -1.361, -1.361) 
(-1.361,1.361,-1.361) 



Notes: (a) van Smaalen et al. (2001) incorrectly reported a (3 + 4)D superspace group P422(<x, fi, y). (b) Wustite, x = 0.098; Yamamoto (1982). (c) Provided for purposes of comparison 
with the other groups with BSG Fm3m; 'published' setting is the symbol from Yamamoto (2005). (dj Ohmasa et al. (1995). (e) Esmaeilzadeh et al. (2001) only discuss the supercentered 
setting. They give two more compositions: Bi 1 _ I Cr I 0 1 5+1 5x with x = 0.05 (a = 0.307) and x = 0.10 (a = 0.2973). (J) Valldor et al. (2000). (g) Withers et al. (1999): Bi^Nb.Oij^ 
(0.06 < x < 0.23). (h) Ling et al. (1998): Bij^Ta^O, 5+I (0.1 < x < 0.25). (i) Esmaeilzadeh et al. (2000). (/) Withers et al. (1990) and Yamamoto (1993). 



q 1 = (0.0641,0.0641,0.151) 
q 2 = (-0.0641, 0.0641, 0.151) 
q 3 = (-0.0641, -0.0641, 0.151) 
q 4 = (0.0641, -0.0641, 0.151). 



q 2 x, 



3_ 

q x, 



4— 1— 

q x = q x 



q 2 x + q 3 x. 



(45) 



This symmetry becomes obvious in the supercentered setting, 
where SSG(3 + d)D shows that there are two symmetry- 

(44) 

equivalent modulation wavevectors, Q 1 and Q 2 , in addition to 
a third wavevector, Q 3 , parallel to the tetragonal axis: 



van Smaalen et al. (2001) incorrectly reported this as a four- 
dimensional modulation (Table 4), but then continued to show 
that the phase transition is accompanied by a lowering of the 
point symmetry and the formation of a multiply twinned 
crystal with a one-dimensional incommensurate modulation in 
each domain. Nevertheless, for the purpose of illustrating a 
fundamental issue of symmetry, we will proceed as though all 
modulation wavevectors would originate in a single domain, 
where the number of symmetry-equivalent modulation 
wavevectors is larger than the dimension of the modulation. 
The modulation in equation (44) is actually three-dimensional, 
because q 4 = q 1 — q 2 + q 3 . Despite this relationship between 
the modulation wavevectors, the tetragonal symmetry requires 
that modulation wavefunctions are symmetric in the four 
arguments 





a 2 * • 



• « 



Figure 3 

Reciprocal plane parallel to (a*, a|) of a hexagonal lattice, (a) Indexing of 
satellite reflections with two modulation wavevectors enclosing an angle 
of 60°. (b) Preferred indexing with modulation wavevectors enclosing an 
angle of 120°. Notice that q' = -q 1 + q 2 . 
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Q 1 = (0.0641,0.0641,0) 

Q 2 = (-0.0641,0.0641,0) 

Q 3 = (0, 0, 0.151). (46) 

The four pairs of satellite reflections as well as the four 
equivalent arguments of the modulation wavefunctions then 
follow as all four equivalent linear combinations of Q 3 with Q 1 
or Q 2 : 

q 1 = Q 1 + Q 3 , q 2 = Q 2 + Q 3 , U1 , 
q 3 = -Ql + Q3, q 4 = _ Q 2 + Q 3_ ^1) 

In accordance with the centering of the superspace lattice, 
diffracted satellite reflections do not appear at ±0^ 
(j = 1, 2, 3), and modulation functions do not contain 
harmonics involving arguments (^x, but only contain linear 
combinations like [Q : x+ Q 3 x] [equation (47)]. 

4.2.2. Superspace symmetry with BSG Fm3m. All known 
compounds with a three-dimensional modulation possess 
cubic symmetry. Wustite, Fe 1 _ x O, is based on an F-centered 
cubic lattice with BSG Fm3m and the simple modulation with 
q 1 = (0.398, 0, 0) (Table 4). The superspace group is 
symmorphic and centerings other than the F-centering, 

(1 , \ , 0, 0, 0, 0) 

(1 , 0, 1 , 0, 0, 0) (48) 
(0,1,1,0,0,0), 

of the BSG do not occur. 

Three compounds have been reported with symmetry 
according to superspace group No. 225.3.215.7 
Fm3m(a, a, a)000(o;, —a, — a)000(— a, a, — a)000. They have 
the same BSG as wustite but different modulation wavevec- 
tors. As with (TaSe 4 ) 2 I (§4.2.1), four symmetry-equivalent 
modulation wavevectors exist. The supercentered setting 
clearly reveals the three-dimensional nature of the modulation 
with q 4 = (-q 1 - q 2 - q 3 ) (Table 4), and 

qi = Q 1 + Q 2 + Q 3 , q 2 = Q 1 _Q 2 _Q 3 , 
q 3 = -Q 1 + Q 2 - Q 3 , q 4 = -Q 1 - Q 2 + Q 3 . 

For Bi 0 8 5Mo 016 O 17 4 the modulation wavevectors of the 
supercentered setting are 

Q 1 = (0.2926, 0, 0) 

Q 2 = (0, 0.2926, 0) (50) 
Q 3 = (0, 0, 0.2926). 

The centering translations of the supercentered setting 
combine the F-center of the basic structure [equation (48)] 
with a so-called 'F-center' among the internal superspace 
coordinates, the latter being defined as 

(0, 0, 0,0,1, Ij 

10, 0, 0, 1 , 0, 1) (51) 
(0, 0,0,1,1,0). 

This can be compared with superspace group No. 225.3.212.5, 
based on modulation wavevectors of the type (0, /S, /3), where 
the supercentered setting again involves modulation wave- 
vectors of the type Q 1 = (/S, 0, 0) [equation (50)], but now 
combines the F-center of the basic structure [equation (48)] 



with an '/-center' among the internal superspace dimensions 
with centering translation (0, 0, 0, \ , \ , £). 

Interestingly, replacing three-valent molybdenum atoms by 
five-valent niobium or tantalum atoms leads to a similar, but 
different structure involving mirror planes with nonzero 
intrinsic translational components along the internal super- 
space dimensions (Table 4). 

4.2.3. Modulation in the /-centered lattice of V 6 Ni 16 Si 7 . 
V 6 Ni 16 Si 7 is a three-dimensionally modulated crystal with 
symmetry based on the cubic /-centered lattice and BSG 
Im3m (Table 4). Withers et al. (1990) report an indexing of the 
electron diffraction based on the modulation wavevectors 

q 1 = 1(1, 1, 4) + e(l, 1, 1) = (0, 0, 1) + a'(l, 1, 1) 

q 2 = 1 (1 , 1, 4) + e(l, 1, 1) = (0, 0, 1) + a'(l , 1, 1) (52) 

q 3 = 1(1, T, 4) + e(T, T, 1) = (0, 0, 1) + a'(T, T, 1) 

where a' = ± + e and q 4 = 1(1, 1, 4) + e(T, 1, 1) = q 1 - q 2 + q 3 . 

Withers et al. (1990) also report the observed reflection 
conditions, but then provide an analysis based on the theory of 
irreducible representations (normal-mode analysis). Yama- 
moto (1993) has assigned to V 6 Ni 16 Si 7 the (3 + 3)D super- 
space group with the tentative symbol Im3m(a, a, 1 + a). 

SSG(3 + d)D shows that such a superspace group does not 
exist. Since symmetry operators are not provided by Withers et 
al. (1990) or Yamamoto (1993), we could not use the f indssg 
tool on SSG(3 + d)D for computing the transformation to the 
standard setting. However, SSG(3 + d)D does show that the 
only possible modulation wavevectors for three-dimensional 
modulations with BSG Im3m are (a, 0, 0), (0, /S, ft) and 
(a, a, a) (Bravais classes 3.208, 3.211 and 3.214, respectively). 
Indeed, the modulation wavevectors can be rewritten as 

q 1 + (1, 1, 0) = (1, 1, 1) + a'(l, 1, 1) = a(l, 1, 1) 

q 2 + (1, 1, 0) = (1, 1, 1) + or*(l, 1, 1) = o(l, 1, 1) (53) 

q 3 + (1, 1, 0) = (1, 1, 1) + a'(\, 1, 1) = a(l, 1, 1), 

where a = l + a' = j + e. Notice that we cannot add the basic 
structure reciprocal vector (0, 0, T) to the modulation wave- 
vectors [equation (19)], because this is a forbidden reciprocal 
vector for the /-centered lattice. Instead, we have added the 
vector (1, 1, 0) to q 1 in order to arrive at a reciprocal vector 
along the diagonal of the cubic unit cell. Of course, this 
goes at the expense of a considerably increased 
length for the modulation wavevectors. Nevertheless, a 
description that respects the symmetry of the problem 
requires these long modulation wavevectors. With the 
new indexing, the non-symmorphic superspace group 
Im3m(a, a, a)q00(a, —a, —a)q0Q(—a, a, — a)000 is obtained, 
which corresponds to No. 229.3.214.8 in SSG(3 + d)D (Table 4). 

5. Incommensurate composite crystals 

Incommensurate composite crystals comprise two or more 
subsystems, each of which has an incommensurately modu- 
lated structure. The basic structures of the subsystems are 
mutually incommensurate, but for all known compounds, any 
pair of subsystems share a common reciprocal-lattice plane of 
their basic structures. The third reciprocal basis vector of one 
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(a) (b) 

Figure 4 

Crystal structure of [Sr] x [TiS 3 ]. (a) Projection showing the two types of 
columns with mutually incommensurate periodicities a 13 for the a 3 axis of 
the first subsystem (TiS 3 ) and a 23 for the a 3 axis of the second subsystem 
(TiS 3 ). (b) Projection along the mutually incommensurate direction 
showing the common basal plane of the hexagonal lattice. Large circles 
denote metal atoms, small circles represent sulfur atoms. Reprinted from 
Figs. 1.4(c) and 1.4(d) in van Smaalen (2007) by permission of Oxford 
University Press (http://www.oup.com). 

subsystem then acts as modulation wavevector for the other 
subsystem, and the other way around. The symmetry of a 
composite crystal is given by a (3 + d)D superspace group, 
while the symmetry of each subsystem is also given by a 
(3 + d)D superspace group. These so-called subsystem 
superspace groups often are different (inequivalent) groups 
according to the definition of superspace group employed in 
de Wolff et al. (1981) and Stokes et al. (2011a). 

The various aspects of the structures and symmetries of 
composite crystals are illustrated by the example of [Sr] x [TiS 3 ] 
(Onoda et al, 1993), where square brackets indicate the 
subsystems. The seemingly non-stoichiometric composition 
with x = 1.132 reflects the incommensurate ratio of the 
volumes of the basic structure unit cells of subsystem 1 (TiS 3 ) 
and subsystem 2 (Sr). [Sr] J [TiS 3 ] is a composite crystal of the 
columnar type, where chains of Sr atoms and columns of TiS 3 
are alternatingly arranged on a two-dimensional hexagonal 
lattice (Fig. 4). The basic structure reciprocal lattices share the 
basis vectors in the basal plane, while the third direction 
(parallel to the chains) is the incommensurate direction: 



a l — a ll — a 21 

n.21 



a 2 — a 12 — a 22 
it 



a 3 — a 13 — q a 4 — a 23 — Q 



(54) 



where a* 3 , for example, denotes the third reciprocal basis 
vector of the first subsystem and q 11 is the first (and in this 
example only) modulation wavevector of the first subsystem 
with q 11 = (a 13 /a 23 K 3 . 

An indexing of all reflections with four integers is obtained 
with the four reciprocal basis vectors M* = [a* , a 2 , a 3 , a^}. 
Along with its modulation wavevector, the reciprocal basis 
vectors of subsystem v (v = 1, 2) are obtained from the four 
reciprocal vectors M* by a (3 + d) x (3 + d) integer matrix W 
(d = 1 in the present example) according to 



3+rf 

k=l 
3+rf 



(55) 



The matrices W extract the basic structure reciprocal basis 
vectors and modulation wavevectors of subsystem v from the 
basis vectors used for indexing. In this sense, W v represents a 
coordinate transformation in superspace between the arbi- 
trarily chosen superspace representation M* and the natural 
subsystem superspace, which is specific to each subsystem. 
Operators of the subsystem superspace group follow as (van 
Smaalen, 1991) 



(56) 



Because reciprocal basis vectors of one subsystem act as 
modulation wavevectors of the other subsystem, W v must be a 
coordinate transformation that mixes the first three dimen- 
sions and the additional dimensions for at least some of the 
subsystems. This coordinate transformation is a forbidden 
transformation when establishing the equivalence of super- 
space groups (Stokes et al, 2011a). Therefore, the subsystem 
superspace groups are generally inequivalent, unless they are 
equivalent by chance, as is the case for the mineral levyclau- 
dite which possesses triclinic symmetry (Evain et al, 2006). 

For [Sr] x [TiS 3 ] equation (54) shows that W 1 is the identity 
matrix. This choice of M* has become a de facto standard for 
composite crystals. It implies a setting where the symmetry of 
[Sr] A .[TiS 3 ] and the symmetry of the first subsystem are 
described by the same superspace group. Onoda et al. (1993) 
give the superspace group P.R3m(0, 0, y)ls, which is found to 
be an alternate symbol for superspace group No. 166.1.22.2, 
R3m(0, 0, y)0s on SSG(3 + d)D. Apart from the ^-centering of 
the hexagonal basic structure unit cell, other centerings in 
superspace do not exist for this lattice. 

Equation (54) leads for the second subsystem to 



W 2 = 



/l 0 0 o\ 

0 10 0 

0 0 0 1 

\0 0 1 0/ 



(57) 



The it-centering of the original setting transforms by W 2 
[equation (56)] into the superspace centering vectors 
(| , | , 0, |) and (| , | , 0, |), which represent an //-type centering 
of the BSG and which has been denoted as the //'-centering of 
the superspace lattice (van Smaalen, 2007). SSG(3 + d)D shows 
that the transformation by W 2 [equation (57)] leads to the 
supercentered setting of the (3 + 1)D superspace group No. 
163.1.23.1, P31c(l,l,y)000. 

The subsystem superspace groups of [Sr] x [TiS 3 ] turn out to 
be inequivalent (3 + 1)D superspace groups, although they are 
of course equivalent as 4D space groups as governed by the 
coordinate transformation W 2 . The case of [Sr] J [TiS 3 ] is 
special as it combines different Bravais lattices of the BSG for 
the subsystems. The rhombohedral lattice with an /^-centering 
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Table 5 

Number of chiral superspace groups in comparison to the number of 
superspace groups. 

Dimension of space or superspace 



Classification 


3 


3 + 1 


3 + 2 


3 + 3 


Bravais classes 


14 


24 


83 


215 


Superspace groups 


230 


775 


3338 


12584 


Chiral superspace groups 


65 


135 


368 


1019 


Fraction that is chiral 


0.283 


0.174 


0.110 


0.081 



Figure 5 

Basal plane of the reciprocal lattice of [Sr] x [TiS 3 ]. A* and are the 
reciprocal basis vectors of the BSG setting of the rhombohedral lattice. 
Filled circles represent Bragg reflections and open circles are the 
positions of extinct Bragg reflections. a*,aj are the reciprocal-lattice 
vectors of the trigonal unit cell in the standard (primitive) BSG setting. 
Extinct Bragg reflections (open circles) only apply to this lattice if the 
supercentered setting (H' setting) is used. 



and the primitive trigonal space group described with an 
//-centered unit cell can be considered as different centerings 
of the hexagonal unit cell, which share a common reciprocal- 
lattice plane perpendicular to the trigonal axis (Fig. 5). 

Other incommensurate composite crystals, including misfit 
layer sulfides (Wiegers, 1996), misfit layer cobalt oxides (Isobe 
et al, 2007) and urea inclusion compounds (van Smaalen & 
Harris, 1996), also exhibit a pairing of two inequivalent 
superspace groups. A detailed analysis of this feature is 
outside the scope of the present overview and will not be 
discussed further here. 

Standard settings and alternate settings of superspace 
groups occur for incommensurate composite crystals by means 
of the same kinds of coordinate transformations as have been 
discussed for modulated crystals. One difference is the 
stronger inclination for employing non-standard settings in the 
case of composite crystals, because the set M* of reciprocal 
basis vectors is preferably chosen to contain reciprocal 
basis vectors of the basic structures of the subsystems. For 
example, for the case of [PbSJj 18 [TiS 2 ] and isostructural 
[Ca 0 85 OH] 1 256 [CoO 2 ] (van Smaalen et al, 1991; Isobe et al, 
2007), M* has been chosen as 



q" 

a 23 



a 2 



a 12 
a 21 



a 22 



(58) 



This results in a mixed setting of the (3 + 1)D superspace 
group as C'2/m(a, 0, 0>0 with C = (§ , \ , 0, £>. SSG(3 + d)D 
shows that this is an alternate setting of superspace group No. 
12.1.7.4 B2/m(0, 0, y)s0. Apart from the trivial transformation 
of the setting of the BSG, the transformation toward the 
standard BSG setting involves the choice of an alternate 
modulation wavevector: 



a 4 = q = a 1 + a 4 = a u + a 21 = a n + q . 



(59) 



The interpretation of a'4 as a reciprocal basis vector of the 
second subsystem is lost in this representation [equation (59)]. 
Therefore, the mixed setting (i.e. not BSG setting or super- 
centered setting) with centering translation (| , \ , 0, |) is 
preferred over the standard setting in the case of these 
composite crystals. 



6. Chiral superspace groups 

Chiral space groups are space groups that may be the 
symmetry of crystals containing chiral molecules. They are of 
particular importance in the life sciences, because all proteins 
and nucleotides are molecules of this type (Lovelace et al, 
2008). 

Chiral space groups are those space groups of which the 
point group contains rotations only (Blow, 2002). Chiral 
superspace groups are then defined as the superspace groups 
for which the three-dimensional point group of the BSG 
contains rotations only (Souvignier, 2003). A list of (3 + d)D 
superspace groups (d = 1, 2, 3) has been generated with this 
criterion and is available on SSG(3 + d)D. It is noticed that the 
fraction of superspace groups that is chiral strongly decreases 
on increasing superspace dimension d (Table 5). We did not 
find a compelling theoretical reason for this feature. But we do 
observe that the number of ways to combine the intrinsic 
translations of the BSG with the intrinsic translations along 
the additional superspace dimensions, or with supercentering 
translations, increases with d. So it appears that the intrinsic 
translations of chiral BSG operations are more restricted in 
the combinations in which they can participate. 

Superspace groups are defined on the basis of equivalence 
relations that only allow coordinate transformations that 
preserve the handedness of the coordinate axes in three- 
dimensional space [det^) = 1; equation (7)], i.e. that 
preserve chirality. This definition leads to pairs of enantio- 
morphic superspace groups in cases where the BSG is an 
enantiomorphic space group, like the (3 + 2)D superspace 
groups No. 76.2.60.2 P^O, 0, ^0(0, 0, y 2 )0 and No. 78.2.60.2 
P4 3 (0, 0, Yi)0(0, 0, y 2 )0. Intrinsic translations along the addi- 
tional superspace dimensions do not give rise to enantio- 
morphic superspace groups. For example, superspace group 
No. 75.2.60.4 P4(0, 0, y^qiO, 0, y 2 )0 is not enantiomorphic 
(q stands for the fractional translation £). Instead of being 
an enantiomorph, P4(0, 0, ^(/(O, 0, y 2 )0 is an alternative 
setting of No. 75.2.60.4, and is transformed into the 
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standard setting by the choice of a different modulation 
wavevector: (0, 0, y[) = (0, 0, — y 1 ) [equation (19)]. 



7. Conclusions 

The computational complexity of finding the transformation 
between two settings of a (3 + <f)D superspace group is 
surprisingly high, especially for d = 2,3. Here an efficient 
algorithm is presented, which either establishes two super- 
space groups to be different superspace groups or determines 
them to be different settings of the same superspace group and 
then provides the transformation between these settings. The 
algorithm has been implemented as an internet-based utility 
called 'findssg', which identifies any user-given (3 + d)D 
superspace group (d = 1, 2, 3) based on the superspace-group 
operators provided, and displays the transformation to the 
standard setting of this superspace group in the SSG(3 + d)D 
tables. 

The algorithm considers coordinate transformations in 
superspace. It is shown that in general such a transformation 
corresponds to one, or a combination, of the following three 
types of transformations in physical space: 

(i) A transformation of the basic structure unit cell. 

(ii) Adding any reciprocal-lattice vector of the basic struc- 
ture to the modulation wavevector [equation (19)]. 

(iii) Replacing originally chosen modulation wavevectors by 
linear combinations of the same [only for d > 2; equation 
(32)]. 

These transformations are illustrated by the analysis of the 
symmetries of a series of compounds with d = 1, 2, 3, 
comparing published and standard settings and discussing 
the transformations between them. It is argued that non- 
standard settings are needed in some cases, while standard 
settings of superspace groups are desirable in other cases. A 
compilation is provided of standard settings of compounds 
with two- and three-dimensional modulations (Tables 2 and 4). 
It appears that several ad hoc notations have been used in the 
literature for (3 + d)Y> superspace groups, especially for d = 2 
and 3. 

For d = 2 superspace groups with trigonal/hexagonal 
symmetry, an angle of 120° between the two modulation 
wavevectors is preferred and is the only correct choice for 
acentric trigonal cases (Table 3). This is the standard setting 
for all relevant Bravais classes in SSG(3 + d)D, in contrast to 
the use of a 60° angle in most published structures (§4.1.3 and 
Table 2). 

The problem of superspace-group settings, including the 
choice of origin, is subtle. Therefore, we strongly advise 
authors to explicitly document for each structure the list of 
symmetry operators (or at least the generators) of the super- 
space group, along with the explicit form of the modulation 
wavevectors as in equation (2). It would also be useful to 
include the number and symbol of the standard setting on 
SSG(3 + d)D for each structure, because this will make it easier 
to check the equivalences of structures and symmetries in 
future studies. 



We thank Gloria Borgstahl for pointing out the importance 
of chiral superspace groups. 
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